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Abstract. Clustering coefficient is an important topological feature of complex 
networks. It is, however, an open question to give out its analytic expression on 
weighted networks yet. Here we applied an extended mean-field approach to investigate 
clustering coefficients in the typical weighted networks proposed by Barrat, Barthelemy 
and Vespignani (BBV networks). We provide analytical solutions of this model and find 
that the local clustering in BBV networks depends on the node degree and strength. 
Our analysis is well in agreement with results of numerical simulations. 
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1. Introduction 

The concept of weight plays an important role in characterizing the features of complex 
systems in nature, since many different correlations in the real world can be represented 
by weighted links. In the past decade, as an emerging network dynamics, weight 
dynamics has attracted a large number of scientific communities [U El El HI El El 
lUISlCin]- It has been shown that many biological, technological, and social systems are 
best described by weighted networks, as the nodes interact with each other with varying 
strength. To describe the variety, a link between two nodes is associated with a weight 
in the weighted network. The weight of a link can represent different properties for 
different networks. For a social network, link weight can indicates the intimacy between 
two individuals [HI [12]; for the Internet, link weight can represent the bandwidth of 
an optical cable between two routers [13]; for the world-wide aviation network, link 
weights can reflect the annual volume of passengers traveling between two airports pQ; 
and for the E. coli metabolic network, link weight can encode the optimal metabolic 
fluxes between two metabolites [H]. 

Previous studies have pointed out the importance of the link weight of networks and 
led to the formulation of a long list of models aimed at incorporating the characteristics 
in weighted networks [U[2l[3l[ll[5l[6l[7l[8l[9l [TQl [15l H6l [13 HH [H] . Among these 
models, the one proposed by Barrat, Barthelemy and Vespignani (BBV) [1], [T9l [20] has 
been deemed to be the pioneer and classic work by most physics communities. The 
main contribution of their work can be ascribed to two aspects: a new strength driven 
preferential attachment mechanism; local weight and strength rearrangements for a new 
link. The results of this model represent many behaviors observed empirically in real- 
world networks, e.g., link weight-degree correlation, node strength-degree correlation, 
and scale-free topology. 

Currently, a primary purpose of the study of complex networks is to understand 
how their dynamical behaviors are influenced by underlying geometrical and topological 
properties [2T1 [22] . Among many fundamental structural characteristics [23], clustering 
coefficient, also known as transitivity, is an important topological feature of complex 
networks, generating considerable attention. In the sociological literature, for example, 
it is referred to as the "network density" . For node i with fcj links, Ci denote the fraction 
of these allowable edges that actually exist. The local clustering coefficient for the whole 
system is given as the average of Ci for each vertex [21]. As is known, average local 
clustering coefficient is one of the two key criterions to judge whether a network or a 
complex system possess the small-world property. Most previous related studies have 
been confined to numerical methods, which is, however, prohibitively difficult for large 
networks because of the limit of time and memory. Hence, there is a need to give out 
a analytical solution of the clustering coefficient toward providing a useful insights into 
topological structures of complex networks. 

Despite the necessity, to the best of our knowledge, the rigorous computation 
for the clustering coefficient of BBV model has not been addressed, contrasting 
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sharply with the successes in other network properties such as degree, weight and 
strength distribution [19]. To fill this gap, in this present paper we investigate this 
interesting quantity analytically. We derive an exact formula for the clustering coefficient 
characterizing the BBV model. The analytic method we adopt in the work, is an 
extended mean field approach. The obtained precise result shows that the average local 
clustering coefficient of BBV network has a power-law decrease with the number of 
nodes. To some extent, our research opens a way to theoretically study the clustering 
coefficient of neutral or nonassortative mixing [25] weighted networks on the one hand. 
On the other hand, our exact solution may give insight different from that afforded by 
the approximate solution of numerical simulation. 



2. Introduction to the model 



The BBV network starts from an initial seed of N vertices fully connected by links 
with assigned weight wq. A new vertex n is added at each time step. This new site 
is connected to m previously existing vertices (i.e., each new vertex will have initially 
exactly m edges, all with equal weight Wq), choosing preferentially sites with large 
strength; i.e., a node i is chosen according to the probability 

n_ = =^-. (i) 

z^ v s v 

where s< is the strength of the node i, defined as Sj = J2j w ij- The weight of each new 
link Lij, is initially set to a given value wq. The creation of this edge will introduce a 
local variations of the edge weight ty i<7 - — > + A Wij , in which node j is one of node 
z's neighbors and node strength of node Sj — > Sj + A^, where this perturbation is 
proportionally distributed among the edges according to their weights. 

= S^. (2) 

Si 

Simultaneously, this rule yields a total strength increase for node i of Si + wq, 
implying that s, — > + Si + w . In this paper, we focus on the case of homogeneous 
coupling with Si = S = const and equivalently set w = 1. After the weights and 
strengths have been updated, the growth process is iterated by introducing a new vertex, 
until the desired size of the network is reached. 

In BBV networks, the dynamical process for s< and fcj can thus be described by the 
following evolution: 

dsi Si(t) Si(t) x Wij(t) 



m 



dki sAt) 

— = m— — — . (4) 

Note that Xljenfi) runs over ^ e neares f neighbors of the node i. Solving these two 
equations with initial conditions fcj(tj) = Sj(tj) = m, one can have 

£ \ 2.S+2 



m ( - ) (5) 
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Figure 1. (Color online) The initial value of the local clustering coefficient Ci(U) for 
(a) 5 — 0.1 and (b) S = 0.2, (averaged over 10 BBV networks). Red circles show the 
analytical solutions given by Eq. (f2"0"|) . Note that the Cj(i) axis is nonlinear. 



P2], and 

Si + 2m5 

ki = -2STT- (6) 

After introducing the BBV model, we now investigate analytically the clustering 
coefficient of all the nodes in the BBV networks. As is known, like BA networks, the 
clustering coefficient in BBV networks rapidly decreases with the network size t. In the 
previous studies, researchers give out numerical calculations to show this mechanism 
while we will offer an exact solution to explain it in what follows. 



3. Analytical expression to clustering coefficients of the model 

Before introducing the way to calculate the clustering coefficient, it is necessary to build 
up a conception, 'what is neutral weighted networks', at first. 

In networks, the degree of a node is defined as the number of links the node has. 
A relevant topological property is the degree-degree correlation, or so-called assortative 
mixing. It is defined as a preference for high-degree vertices to attach to other high- 
degree vertices and vice versa [25l [26] . By contrast, disassortative mixing is defined as a 
preference for high degree vertices preferentially connect with low degree ones and vice 
versa. However, a network with neutral degree correlation means neither assortative nor 
disassortative. A random network, where links are randomly connected between nodes, 
is an example of neutral networks. 

Along with the degree correlations, we introduce the weighted degree correla- 
tions pQ, defined as 

1 N 

^nn,i = 7" ^ ClijWijkj. (7) 
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Figure 2. (Color online) The local clustering coefficient Ci(t) as a function ti for 
(a) S = 0.1 and (b) 5 — 0.2. Results are obtained by numerical simulations on BBV 
networks (averaged over 10 networks). Peachblow lines show the analytical solutions 
given by equation (|19p. and black curves are the smoothed data of simulation data 
(gray circles). 



The k™ ni thus measures the effective affinity to connect with high- or low-degree 
neighbors according to the magnitude of the actual interactions. As well, the behavior of 
the function k™ n (k) marks the weighted assortative or disassortative properties. When 
k™ n (k) increases with k, it means that nodes have a tendency to connect to nodes with 
a similar or larger degree. In this case the network is defined as assortative. In contrast, 
if k™ n (k) is decreasing with k, which implies that nodes of large degree are likely to 
have near neighbors with small degree, then the network is said to be disassortative. 
Naturally, if correlations are absent, it can be ascribed to weighted neutral networks. 

Of particular interest is the BBV model mentioned in section [2] can give rise to 
neutral networks, for small S [19]. However, for increasing 5, the disassortative character 
and a power-law behavior of k nn (k) emerge. Fortunately, their weighted assortativity 
k™ n (k) is a constant, that is, they belong to weighted neutral networks. In this case, we 
redefine the linking probability of node i and j, 

Pij = = 2m(i + <y)twv (8) 

for the sake of extending the mean field method to this particular network for large 5. 

In a BBV network, consider consequently a certain node i at time step t. Note that 
we only consider the case m > 1 in that the BBV network is close to tree-like networks 
when m = l. As is known to us all, for a tree network, the clustering coefficient in the 
network is equal to zero. By the definition [21], the clustering coefficient Ci depends on 
two variables and fcj. Its definition is Ci = k .^ i _ 1 ^ , where E^ is the number of links 
between the k t neighbors of i. 
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Since in the BBV model only new nodes may create links, the coefficient Cj changes 
only when its degree k{ changes, i.e., when new nodes create connections to i and 
x G (0, m — 1) of its nearest neighbors. The appropriate equation for changes of Cj is 
then 

-^ i = X)PteA Gte , (9) 

x=0 

where Ac ix denotes the change of the clustering coefficient when a new node connects 
to the node i and to x of its first neighbors as well, whereas pi X describes the probability 
of this event. 

Notice that, in the equation (jUJ), Ac lx represents the difference between clustering 
coefficients of the same node i calculated after and before a new node attachment 

2(Ei + x) 2Ei 2x 2d 

Ci * ~ kjkt + 1) " ki(ki - 1) ~ ki{ki + 1) ~ k+l ' ( } 

where A Cix > 0. Also, the equation (|T0|) gives out the restricted region of our method. 
The probability pi X is a product of two factors. 

p ix = P-(l - P)— i-. (11) 

The first factor is the probability of a new link to end up in z, which is given by 
equation ([1]). The second one is the probability that x among the rest of (m — 1) new 
links connect to the nearest neighbors of i. It is equivalent to the probability that (m— 1) 
Bernoulli trials with the probability P for x successes. So that, 



E w s« 2m(l + «J)*" 
Replacing the sum Ejen(i) by an integral, one obtains 



(12) 



ien(i) Jl 



where, 



in which Wij is the weight of the link Lij. In BBV networks, 

8 

Wij(t) = i^-j , (15) 
where t^ = max(i, j). Hence, using the equation (fT3l) and one can find the probability 

4,5-1 / 66+1 23+1 26+1 

mt 2S + 2 / „4i+l i 25+2 f 2S+2i 2S+2 

i 25+2 . (16) 



4(1 + 8)5 
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Now, inserting equation flSD, flSD, (fTUj) and (TTT|) into equation 
dCi m(25 + l)t^d 



one can obtain 



dt 



2.5+1 



26 + 1 



-3 
£2,5+2 



1 + 5)[mt^ + (2m5 + 25 + l)tf 6+2 '\ 
m(m - 1)(2S+ 1) 



;i7) 



In the derivation process, we simplify 



26+1 



m t2s+2 + (2m<5 + 26 + l)i 

26+1 26+1 

mt^+2 + 2m5i 2S + 2 



2S + 1 
25+2 



to make the analytical integral easily executed. Solving the equation for Cj one gets 



Ci{t) 



Bi + 



2S+1 



mt2s+2 + (2m5 + 25 + l)t 



25+1 
2,5 + 2 



1 -2 



3d 



3,5 



m 



4(l + <5)5 2 



2t^t i 6+1 - 5+1 - 6t^+r 

6 



(19) 



where Bi is an integration constant for the node i and determined by the initial condition 
Ci(ti) that describes the clustering coefficient of the node % exactly at the moment of 
its attachment t{. Note that Ac ix > in equation (TTUl) . we have p > Cj. Thus, in 
this article, it is worth noticing that 5 has to be confined in the region [0, 1/4). The 
restriction is given by our approximation in equation ( JTHl) . We believe this drawback 
will be addressed in future work. 

To clarify the calculation of Cj(tj), we show another form of the definition of 

Ci [22J Elj, 

number of triangles connected to vertex i 



It can be rewritten as 



number of triples centered on vertex i 



Ci{ti 



E ? E v Pij Viv Pjt 



m 



1-t 



-6 
5+1 



1-t 



25 
5+1 



5-1 
>.5+l 



(20) 



2(7) 16(m- l)(5+ 1)S 2 

Figure Q] shows the prediction of the equation (1201) in comparison with numerical 
results. For small values of £j the numerical data differ from the theory significantly. 
This incongruity is caused by the formula for the probability of a connection pij in 
equation (jSj) used three times in equation (12"0"|) . which holds only in the asymptotic 
region t { — > oo. 

Taking into account the initial condition Cj(ij) from equation (1201) . one can obtain 
Bi of a given node i 



m\2mb + 26 + m + l) 2 1 - t 



B, 



1 - t 



-26 
5+1 



33 
t/ +1 



16(m- !)(<$ + 1)5 2 



rrr 



2a 
«+i 



3 

7t/ +1 



4(l + <5)<5 2 



(21) 
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Figure 3. (Color online) The clustering coefficient C of a whole BBV network as a 
function of the network size t for different values of the 8. Results are obtained by 
numerical simulations on BBV networks (averaged over 10 networks). The solid lines 
are numerical solutions given by equation (I22p . 



To obtain the average clustering coefficient C of the whole network, the expression 
equation (j2ip has to be averaged over all nodes within a network. Admittedly, it is 
difficult to find an exact analytic result of C uniformly but we can give out numerical 
solutions of BBV networks by the definition 

C = 22^1 . (22) 

In the figure ([3]), one can easily find out that both the numerical and theoretical results 
exhibit a power-law decay as the size of networks growing. Hence, when t tends to 
infinity, C is getting close to zero. As shown in the figure, the power-law exponent 
depends on the the tunable parameter 5. For large 5, the decreasing velocity is relatively 
slow. 
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4. Conclusion 

To sum up, the clustering coefficient in complex systems is an important topological 
feature of complex systems. It has a profound impact on a variety of crucial fields, such 
as epidemic processes, network resilience and finding community structure, to name 
but a few. In this paper, we have derived analytically the solution for the clustering 
coefficient of the Barrat et al.'s model [19] which has been attracting much research 
interest. We found that the average clustering coefficient decays with the size of networks 
growing, following a power-law. Our analytical technique could guide and shed light 
on related studies for weighted networks by providing a paradigm for calculating the 
clustering coefficient. 
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